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Abstract 

Let G(n, M) be the uniform random graph with n vertices and M edges. Erdos and Renyi 
(1960) conjectured that the hmiting probabihty 

hm Pr{G{n, ?) is planar} 

n— ^oo 

exists and is a constant strictly between and 1. Luczak, Pittel and Wierman (1994) proved 
this conjecture and Janson, Luczak, Knuth and Pittel (1993) gave lower and upper bounds 
for this probability. In this paper we determine the exact Hmiting probability of a random 
graph being planar near the critical point M — n/2. For each A, we find an exact analytic 
expression for 

p(A) = lim Pr|G(^7i, f(l + An~^/^)) is planar | . 

In particular, we obtain p{0) « 0.99780. We extend these results to classes of graphs closed 
under taking minors. As an example, we show that the probability of G{n, ^) being series- 
parallel converges to 0.98003. For the sake of completeness and exposition we reprove in a 
concise way several basic properties we need of a random graph near the critical point. 



We dedicate this paper to the memory of Philippe Flajolet. 

1 Introduction 

The random graph model G(n, M) assigns uniform probabihty to graphs on n labehed vertices 
with M edges. A fundamental result of Erdos and Renyi [7] is that the random graph G{n, M) 
undergoes an abrupt change when M is around n/2, the value for which the average vertex degree 
is equal to one. When M — cn/2 and c < 1, almost surely the connected components are all of 
order O(logn), and are either trees or unicyclic graphs. When M — cn/2 and c > 1, almost surely 
there is a unique giant component of size 8(n). We direct to reader to the reference texts |3] and 
[Hj for a detailed discussion of these facts. 

We concentrate on the so-called critical window M = f (1 + An^^/^^), where A is a real number, 
identified by the work of Bollobas [21 [3] . Let us recall that the excess of a connected graph is the 
number of edges minus the number of vertices. A connected graph is complex if it has positive 
excess. As A — >■ — oo, complex components disappear and only trees and unicyclic components 
survive, and as A ^ -l-oo, components with unbounded excess appear. A thorough analysis of 
the random graph in the critical window can be found in [13' and [16' , which constitute our basic 
references. 
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For each fixed A, we denote the random graph G (n, f (1 + An^^/^)) by G(A). The core C(A) of 
G{X) is obtained by repeatedly removing all vertices of degree one from G(A). The kernel K{X) is 
obtained from C(A) by replacing all maximal paths of vertices of degree two by single edges. The 
graph G(A) satisfies almost surely several fundamental properties, that were established in [IS] by 
a subtle simultaneous analysis of the G{n, M) and the G(n,p) models. 

1. The number of complex components is bounded. 

2. Each complex component has size of order n^/^, and the largest suspended tree in each 
complex component has size of order n^/^. 

3. C{X) has size of order n^/"^ and maximum degree three, and the distance between two vertices 
of degree three in C(A) is of order n^/'^. 

4. K{X) is a cubic (3-regular) multigraph of bounded size. 

The key property for us is the last one. It implies that almost surely the components of G(A) are 
trees, unicyclic graphs, and those obtained from a cubic multigraph K by attaching rooted trees 
to the vertices of K, and attaching ordered sequences of rooted trees to the edges of K. Some 
care is needed here, since the resulting graph may not be simple, but asymptotically this can be 
accounted for. 

It is clear that G(A) is planar if and only if the kernel K{X) is planar. Then by counting 
planar cubic multigraphs it is possible to estimate the probability that G(A) is planar. To this 
end we use generating functions. The trees attached to K{X) are encoded by the generating 
function T{z) of rooted trees, and complex analytic methods are used to estimated the coefficients 
of the corresponding series. This allows us to determine the exact probability 

p(A)= lim Pr|G(n, ;^(l + An-i/3)) is planar) . 

In particular, we obtain p(0) ~ 0.99780. 

This approach was initiated in the seminal paper by Flajolet, Knuth and Pittel |H|, where the 
authors determined the threshold for the appearance of the first cycles in G(n, M). A basic feature 
in [5j is to estimate coefficients of large powers of generating functions using Cauchy integrals and 
the saddle point method. This path was followed by Janson, Knuth, Luczak and Pittel [13], 
obtaining a wealth of results on G(A). Of particular importance for us is the determination 
in [13] of the limiting probability that G(A) has given excess. The approach by Luczak, Pittel 
and Wierman in [IB] is more probabilistic and has as starting point the classical estimates by 
Wright [T^ on the number of connected graphs with fixed excess. The range of these estimates 
was extended by BoUobas [2] and more recently the analysis was refined by Flajolet, Salvy and 
Schaeffer [9] , by giving complete asymptotic expansions in terms of the Airy function. 

The paper is organized as follows. In Section [2] we present the basic lemmas needed in the 
sequel. Except for the proof of Lemma [T] the paper is self-contained. Lemmas [2] to [5] are proved 
in [13] with a different presentation; for the sake of completeness and exposition we provide 
shorter and hopefully more accessible proofs. In Section [3] we compute the number of cubic 
planar multigraphs, suitably weighted, where we follow [151 . In Section |4] we compute the exact 
probability that the random graph G(A) is planar as a function of A. We generalize this result by 
determining the probability that G(A) belongs to a minor-closed class of graphs in several cases 
of interest. 

We close this introduction with a remark. The problem of 2-satisfiability presents a striking 
analogy with the random graph process. Given n Boolean variables and a conjunctive formula 
of M clauses, each involving two literals, the problem is to determine the probability that the 
formula is satisfiable when M grows with n. The threshold has been established at Af = n and 
the critical window is also of width v?/^] see [S]. However the exact probability of satisfiability 
when the number of clauses is n{l -t- An~^/^) has not been determined, and appears to be a more 
difficult problem. 
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2 Preliminaries 



All graphs in this paper are labelled. The size of a graph is its number of vertices. A multigraph 
is a graph with loops and multiple edges allowed. 

We recall that the generating function T{z) of rooted trees satisfies 

T(z) =ze^(^\ 

Using Lagrange's inversion [TU], one recovers the classical formula n"~^ for the number of rooted 
labelled trees. The generating function for unrooted trees is 



U{z) = T(z) - 




This can be proved by integrating the relation T{z) = zU'{z), or more combinatorially using the 
dissimilarity theorem for trees |17] . 

A graph is unicyclic if it is connected and has a unique cycle. Unicyclic graphs can be seen 
as an undirected cycle of length at least three to which we attach a sequence of rooted trees. 
Since the directed cycle construction corresponds algebraically to log(l/(l — T{z)) (see [TOj), the 
generating function is 

Viz)^Ulo,^^-Tiz)-^). 
^ ' 2\ ^ l-r(z) ^ ' 2 J 

Graphs all whose components are unicyclic are given by the exponential formula: 

yr^T(i) 

The following result, which is fundamental for us, is proved in |16l Theorem 4] by a careful 
analysis of the structure of complex components in G(A). We say that a property V holds asymp- 
totically almost surely (a.a.s.) in G{n,M) if the probability that V is satisfied tends to one as 
n-^oo. Recall that G(A) = G (n, f (1 + An-^/^)). 

Lemma 1. For each X, the kernel of G{X) is a.a.s. a cubic multigraph. 

Given a cubic multigraph M with a loops, b double edges and c triple edges, define its weight as 

w{M) = 2-''2-''6~^ 

This weight (called the compensation factor in [13]), has the following explanation. When we 
substitute edges of the kernel by sequences of rooted trees, a loop has two possible orientations 
that give raise to the same graph. A double (triple) edge can be permuted in two (six) ways, 
again producing the same graph. From now on, all multigraphs we consider are weighted, so that 
we omit the qualifier. The following lemma is proved in |13j using a combination of guessing and 
recurrence relations. The next proof is already contained in [4l Chap. 2]. 

Lemma 2. The number Er of cubic multigraphs with 2r vertices is equal to 

(6r)! 
(3r)!23'-62'-' 

Proof. A cubic multigraph can be modeled as a pairing of darts (half-edges), 3 for each vertex, 
with a total of 6r darts. The number of such pairings is (6r)!/((3r)!23''). However, we have to 
divide by the number G^*" of ways of permuting each of the 2r triples of darts. The weight takes 
care exactly of the number of times a cubic multigraph is produced in this process. □ 

The next result is essentially proved in |13j using several algebraic manipulations. Here we 
present a concise proof. We denote by [z"]A(z) the coefficient of z" in the power series A{z). 
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Lemma 3. The number g(n, M, r) of simple graphs with n vertices, M edges and cubic kernel of 
size 2r satisfies 

and 

g{n, M, r) > n! [z"] ^i'^'''^'^ ^vi^) ^^^]'[,, . 

' ' '- ^ \n-M + r)\ (2r)! (1-T(z))3'- 

Proof. Such a graph is the union of a set of s unrooted trees, a set of unicychc graphs, and a cubic 
multigraph K with a rooted tree attached to each vertex of K and a sequence (possibly empty) 
of rooted trees attached to each edge of K. Let us see first that s — n — M + r . Indeed, the final 
excess of edges over vertices must be M — n. Each tree component contributes with excess —1, 
each unicyclic component with excess 0, and K (together with the attached trees) with excess r. 
Hence M — n — —s + r. 

The first two factors J7(z)"^^^+'7 (n— Af+r)! and e^'-^-' on the right-hand side of the inequalities 
encode the set of trees and unicyclic components. The last part encodes the kernel K . It has 2r 
vertices and is labelled, hence the factor £'r/(2r)!; the weighting guarantees that each graph 
contributing towards g{n, M, r) is counted exactly once. The trees attached to the 2r vertices give 
a factor T{z)^^. The sequences of trees attached to the 3r edges give each a factor 1/(1 — T{z)). 
However, this allows for the empty sequence and the resulting graph may not be simple, so we get 
only an upper bound. To guarantee that the final graph is simple we take sequences of length at 
least two, encoded by r(z)^/(l — T{z)) (length one is enough for multiple edges of K, but length 
two is needed for loops). Since this misses some graphs, we get a lower bound. □ 

The following technical result is essentially Lemma 3 from |13l. We reprove it here for com- 
pleteness in a simplified version tailored to our needs (see also the proof of Theorem 5 in jS,). 

Lemma 4. Let M = ^(1 + Xn^^/^). Then for any fixed a and integer r > Q we have 

^iXjf \n-M + r)\{l-T{z)fr^ - V2^A (3r + , , Aj yi + U^^^^^^jj (1) 
uniformly for |A| < n^/^^, where 

' ^) = 5feTI)75 /fc!r((y+l~2fc)/3) ■ 

Proof. The proof is based on relating the left-hand side of Equation ([ij to the integral represen- 
tation of ^(y,A) defined in [13 Equation(10.7)]: 

2711 Jyi 

where K{X, s) is the polynomial 

(^+Am-A)_.3^,2 ^3 

and H is a path in the complex plane that consists of the following three straight line segments: 

( -e-"/3 for - 00 < t < -2, 
s{t) ^ I 1 -|-itsin7r/3, for - 2 < t < +2, 
[ e+^/^t, for + 2<t<+oo. 

The constant term in the left-hand side of ([T]) is estimated using Stirling's formula, getting 

^'-^e-^'/^-^f^'-(l + 0(^]). (3) 



n! 1 ^ — 2"-*^+'' _,3, 



^(^)^ {n-M + r)! 



4 



The coefBcient of [2"] in Equation ([!]) is estimated by means of a contour integral around z = 0, 
using the expressions of U{z) and V{z) in terms of T{z): 

1 I frj^f . Tizfy-^'^'' T{zYe-n-)l^-n^)"l^ dz 
1 (z) 



27ri 7 V 2 y (1 - r(z))3'^+i/2 ■ 

We make the change of variable u = T{z), whose inverse is z = we"", and we obtain 

?^/.W.-'<"'^^, (4) 
where the integrand is split into a smooth function 

9(-) = (l„,)3.-l/2 

and a large power involving 

f M\ 1 
h{u) = It — 1 — logu — 1 I log 



The contour path in Equation Q should be such that |m| < 1. As remarked in [T3] (see also [5]), 
the function h{u) satisfies h{l) — h'{l) = 0. Moreover, precisely at the critical value M = n/2 we 
also have h"{l) — 0. This triple zero shows up in the procedure used in [13^ when estimating Q 
for large n by means of the saddle-point method. Notice that h{u) is singular at u = 1, due to the 
singularity at z — of T{z). 

Let V = n^^^^, and let a be the positive solution to 

A — aT^ — a. 

This choice is necessary in order to get precise bounds for the tail estimates that appear using 
the saddle-point method. Following the proof of [131 Lemma 3], we evaluate Q on the path 
u = e~("+**)'', where t runs from —nn^/^ to nn^/'^. That is, 

The main contribution to the value of this integral comes from the vicinity of i = 0. The magnitude 
of e'**^"^ depends on the real part ^h{u). Observe that 3fJ/i(e~'-"+**-''^) decreases as \t\ increases, 
and that |e"''("^| has its maximum on the circle u — e~("+**)'' when t — 0. We write s = a + it. 
Analyzing n/i(e~*'^), we have 

nh{e-n = ^ + ^ + 0{iX's' + s^),y), 
uniformly in any region such that \siy\ < log 2. For the function g{u)^ we have 



If f{u) = g(u)e"'''^"'' is the integrand of (|4j), we have 
when s = 0{n^^^^). Finally, 

^ ' /(m)— = e-3/4^3/2-3r^(3^ + i. A) + 0(j.5/2-3rg-A-V6^3r./2+l/4)^ 



27ri / u 



where the error term has been derived from those already in |13j . The proof of the lemma is 
completed by multiplying ([3| and Q, and canceling equal terms. □ 
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It is important to notice that in the previous lemma the final asymptotic estimate does not 
depend on the choice of a. The next result is a direct consequence and can be found as Formula 
(13.17) in 13] . 

Lemma 5. The limiting probability that the random graph G(A) has a cubic kernel of size 2r is 
equal to 

V2^erA{3r+ i, A), 

where = £'f./(2r)! and A{y, A) is as in the previous lemma. 
In particular, for A = the limiting probability is 

Ms [3 J ^'■(2r)!- 

Proof. Using the notation of Lemma |3j the probability for a given n is by definition 



g{n, M,r) 




Lemma [3] gives upper and lower bounds for this probability, and Using Lemma [4] we see that both 
bounds agree in the limit and are equal to 

^V2^A(3r+i,A), 

thus proving the result. A key point is that the discrepancy between the factors T(z)^'' and T(z)^'' 
in the bounds for g{n, M, r) does not affect the limiting value of the probability. □ 

Notice that if we replace the Cr by the numbers gr arising by counting planar cubic multigraphs, 
we obtain immediately the probability that G(A) has a cubic planar kernel of size 2r. Since G(A) 
is planar if and only if its kernel is planar, we can use this fact to compute the probability of G(A) 
being planar. But first we must compute gr. 



3 Planar cubic multigraphs 

In this section we compute the numbers Gr of cubic planar multigraphs of size 2r. The associated 
generating function has been obtained recently by Kang and Luczak |15) (generalizing the enu- 
meration of simple cubic graphs in [I]), but their derivation contains some minor errors. They do 
not affect the correctness of [15], since the asymptotic estimates needed by the authors are still 
valid. However, for the computations that follow we need the exact values. The next result is 
from 15], with the corrections mentioned below. All multigraphs are weighted as in the previous 
section. 

Lemma 6. LetGi{z) be the generating function of connected cubic planar multigraphs. ThenGi{z) 
is determined by the following system of equations: 

Sz"^ = D{z)+C{z) 



dz 



z^ 



Biz) - -iD{z) + Ciz)) + - 

C{z) = S{z) + P{z) + H{z) + B{z) 



z^ 

S{z) ^ C{zf ~ C[z)S{z) 

P{z) = ,2c^,) + ^,^cizr + ^ 

2(1 + C{z))H{z) = u{z){l-2u{z)) -u{z){l-u{z)f 

z^{C{z) + lf = u{z){l-u{z)f. 
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The generating functions B{z), C{z), D{z), S{z), P{z) and H{z) correspond to distinct fami- 
lies of edge-rooted cubic planar graphs, and u{z) is an algebraic function related to the enumeration 
of 3-connected cubic planar graphs (dually, 3-connected triangulations) . 

The corrections with respect to [13] are the following. In the first equation a term — 7z^/24 
has been removed. In the second and sixth equations we have replaced a term by In 
the fourth equation we have removed a term — For the combinatorial interpretation of the 
various generating functions and the proof of the former equations we refer to |15j . Notice that 
eliminating u{z) from the last two equations we obtain a relation between C{z) and H{z). This 
relation can be used to obtain a single equation satisfied by C(z), eliminating from the remaining 
equations. We reproduce it here in case the reader wishes to check our computations. 

1048576 z^ -f 1034496 z^ - 55296 z^ + (9437184 z^ + 6731264 z^ - 1677312 z^ -I- 55296) C+ 

(37748736 z^ + 18925312 - 7913472 z^ + 470016) C^-H 

(88080384 4-30127104 - 16687104 z^ + 1622016) 

(132120576 z^ -|- 29935360 z^ - 19138560 z^ -f 2928640) 

(132120576 z^ -|- 19314176 z^ - 12429312 z^ + 2981888) 

(88080384 z^ + 8112384 z" - 4300800 z^ + 1720320) 

(37748736 z^ + 2097152 z" - 614400 z^ + 524288) C'^+ 

(9437184 z^ + 262144 z^ + 65536) + 1048576 C^z^ = 0. 



The first terms are 



^, , 2 25 4 59 6 11339 c 



This allows us to compute -B(z), D{z), S(z), P{z) and H{z), hence also Gi(z). The first coefficients 
of Gi(z) are as follows. 

^ , , 5 2 5 4 121 1591 o 

Gi z) = — z^ H z^ H z^ H z^ H 

^ ' 24 16 128 384 

Using the set construction, the generating function G(z) for cubic planar multigraphs is then 

cuz) z^"' . 5 2 385 4 83933 g 35002561 o 

G(z) = e'^i^^'' = > Gr-, — TT = iH -\ z^H z^ + z^ -\ , (5) 

^ (2r)! 24 1152 82944 7962624 ' ^' 

r=0 ^ ' 

where Gr is the number of planar cubic multigraphs with 2r vertices. This coincides with the gen- 
erating function for all cubic (non- necessarily planar) multigraphs up to the coefficient of z*. The 
first discrepancy is in the coefficient of z^. The difference between the coefficients is 1/72 = 10/6!, 
corresponding to the 10 possible ways of labelling if 3, 3, the unique non-planar cubic multigraph 
on six vertices. 



4 Probability of planarity and generalizations 

Let G be a graph with a cubic kernel K . Then clearly G is planar if and only if K is planar, 
and we can compute the probability that G(n, M) is planar by counting over all possible planar 
kernels. 

Theorem 7. Let gr(2ry. be the number of cubic planar multigraphs with 2r vertices. Then the 
limiting probability that the random graph G{n, M = §(1 + Xn^^^^)) is planar is 



p{\)^Y.'^9rA{ir+\,\). 



r>0 



In particular, the limiting probability that G{n, f ) is planar is 

/2/4y ^ 

3 [3) ^'■(27^ 



^(0)-E\/If|)'^5.7^« 0.99780. 



r>0 
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Proof. The same analysis as in Section [2| shows that ^/Qjt grA[?>r + 5, A) is the probabihty that the 
kernel is planar and has 2r vertices. Summing over all possible r, we get the desired result. □ 

As already mentioned, Erdos and Renyi |7j conjectured that p(0) exists and < p(0) < 1. This 
was proved in [IS], showing that p{\) exists for all A and that < p(A) < 1. The bounds in [T3] 
for p(0) are 

0.98707 < p(0) < 0.99977, 

obtained by considering connected cubic multigraphs with at most six vertices. We remark that 
Stepanov ll8j showed that p{X) < 1 for A < (without actually establishing the existence of the 
limiting probability). The function p(A) is plotted in Figure 1. As expected, p{X) is close to 1 when 
A — > —00 and close to when A — >■ 00. For instance, p(— 3) « 1 — 1.02 • 10"'' and p(5) « 4.9 • 10"^. 

Besides planar graphs, one can consider other classes of graphs. Let Q he a, class of graphs 
closed under taking minors, that is, if is a minor of G and G E Q, then H E Q. If iJi, • • • , Hk 
are the excluded minors of Q, then we write Q = Ex{Hi, . . . , Hf^). (By the celebrated theorem of 
Robertson and Seymour, the number of excluded minors is finite, but we do not need this deep 
result here). The following result generalizes the previous theorem. 

Theorem 8. Let Q = Ex(Hi, . . . , Hk) and assume all the Hi are 3-connected. Let hr{2ry. be the 
number of cubic multigraphs in Q with 2r vertices. Then the limiting probability that the random 
graph G{n,M = f (1 + Xn~^/^)) is in g is 

pg(\) = ^\/2^/i^A(3r+ 5, A). 

r>0 

In particular, the limiting probability that G{n, ^) is in Q is 
Moreover, for each A we have 

<pg(A) < 1. 

Proof. If all the Hi are 3-connected, then clearly a graph is in Q if and only its kernel is in Q. The 
probability pg{\) is then computed as in Theorem [t] It is positive since Q contains all trees and 
unicyclic graphs, which contribute with positive probability (although tending to as A — ?> 00). 
To prove that it is less than one, let t be the largest size of the excluded minors Hi. By splitting 
vertices it is easy to construct cubic graphs containing Kt+i as a minor, hence G{X) contains Kt+i 
as a minor with positive probability (alternatively, see the argument at the end of [IFj). It follows 
that 1 -pg(A) > 0. □ 

In some cases of interest we are able to compute the numbers h,. explicitly. Let Q = Ex(/'ir4) be 
the class of series-parallel graphs. The same system of equations as in Lemma |6] holds for series- 
parallel graphs with the difference that now H{z) — (this is due to the fact that there are no 
3-connected series-parallel graphs). The generating function for cubic series-parallel multigraphs 
can be computed as 

^ , , 5 2 337 4 55565 . 15517345 . 

GsJz) = IH H z^ H z^ -I z^ H 

' 24 1152 82944 7962624 

For instance, [z^](G'(2;) — Gsp{z)) = corresponding to the fact that K4 is the only cubic 
multigraph with 4 vertices which is not series-parallel. The limiting probability that G(n, ^) is 
series-parallel is 

Psp(O) ~ 0.98003. 

See Figure 1 for a plot of Psp(A). 
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Figure 1: The probability of G(A) being planar and of being series-parallel are both plotted for 
A € [— 1..4]. The function on top corresponds to the planar case. 



The class Ex(i4'4, i^2,3) of outerplanar graphs does not fall directly under this scheme, since 
K2.3 is not 3-connected, but adapting the equations in Lemma [6] (in particular the parallel de- 
composition encoded by P{z)) it is possible to enumerate exactly cubic outerplanar multigraphs. 
The first terms in the generating function are 

^ , , 5 2 337 4 55565 . 14853793 o 

Gout(^) = IH 2^ H H -\ z^ -{ 

^ ' 24 1152 82944 7962624 

The first discrepancy with Gsp{z) is at z*, corresponding to the graph K2^3 with either a loop or 
a double edge attached at the vertices of degree two. The probability of being outerplanar is 

P out (0) « 0.97979. 

We do not plot Pout(A) in Figure 1 since it is too close to Psp(A) to see a clear distinction. 

As another example, consider excluding K^ ^. Since the only 3-connected non-planar graph 
in Ex(i^^3 3) is if 5, which is not cubic, the limiting probability of being in this class is exactly 
the same as of being planar, although Ex(iir3_3) is exponentially larger than the class of planar 
graphs P3]. But excluding the graph if ^3, obtained by adding one edge to if 3, 3, does increase the 
probability, since if3,3 is in the class and is cubic and non-planar (the probability is computable 
since the 3-connected graphs in Ex(if;^3) are known [H]). Other classes such as Ex(if5 — e) or 
Ex(if3 X if2) can be analyzed too using the results from [T^ . 

It would be interesting to compute the probability that G{\) has genus g. For this we need to 
count cubic multigraphs of genus g (orientable or not). We only know how to do this for g — Q, 
the reason being that a 3-connected planar graph has a unique embedding in the sphere. This is 
not at all true in positive genus. It is true though that almost all 3-connected graphs of genus g 
have a unique embedding in the surface of genus g (see [Sj). This could be the starting point for 
the enumeration, by counting first 3-connected maps of genus g (a map is a graph equipped with 
a 2-cell embedding). But this is not enough here, since we need the exact numbers of graphs. 
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